Magnetic moment manipulation by a Josephson current 
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We consider a Josephson junction where the weak-link is formed by a non-centrosymmetric fer- 
romagnet. In such a junction, the superconducting current acts as a direct driving force on the 
magnetic moment. We show that the a.c. Josephson effect generates a magnetic precession provid- 
ing then a feedback to the current. Magnetic dynamics result in several anomalies of current-phase 
relations (second harmonic, dissipative current) which are strongly enhanced near the ferromagnetic 
resonance frequency. 



Many interesting phenomena have been observed re- 
cently in the field of spintronics: the spin-dependent elec- 
tric current and inversely the current-dependent magne- 
tization orientation (see for example [l|, Moreover, 
it is well known that spin-orbit interaction may be of 
primary importance for spintronic, namely for systems 
using a two-dimensional electron gas [s]. In the super- 
conductor/ferromagnet /superconductor (S/F/S) Joseph- 
son junctions, the spin-orbit interaction in a ferromagnet 
without inversion symmetry provides a mechanism for a 
direct (linear) coupling between the magnetic moment 
and the superconducting current Q. Similar anomalous 
properties have been predicted for Josephson junctions 
with spin-polarized quantum point contact in a two di- 
mensional electron gas [5]. S/F/S junctions are known 
to reveal a transition to 7r-phase, where the supercon- 
ducting phase difference Lp in the ground state is equal 
to TT [6|]. However, the current-phase relation (CPR) in 
such a TT-junction has a usual sinusoidal form I — Ic sin (p, 
where the critical current Ic depends in a damped oscilla- 
tory manner on the modulus of the ferromagnet exchange 
field. In a non-centrosymmetric ferromagnetic junction, 
called hereafter (/jq -junction, the time reversal symme- 
try is broken and the CPR becomes I — Ic sin {ip — ipo), 
where the phase shift (fo is proportional to the magnetic 
moment perpendicular to the gradient of the asymmetric 
spin-orbit potential [4]. Therefore, manipulation of the 
internal magnetic moment can be achieved via the super- 
conducting phase difference {i.e. by Josephson current). 

In the present work we study theoretically the spin dy- 
namics associated with such (/3o-junctions. Though there 
is a lot of experimental progress in studying the static 
properties of S/F/S junctions, little is known about the 
spin-dynamics in S/F systems. Note here the pioneering 
work [71| where a sharpening of the ferromagnetic reso- 
nance was observed below the superconducting transition 
in Nb/Ni8oFe2o system. Theoretically, the single spin dy- 
namics interplay with a Josephson effect has been studied 
in [1, [13, El ■ More recently, the dynamically induced 
triplet proximity effect in S/F/S junctions was studied in 
ITsj. while the junctions with composite regions (in- 
cluding several F regions with different magnetization) 
were discussed in [ij, ll5|. Here we consider a simple 
S/F/S (^Q-junction in a low frequency regime hujj <C Tc 



(uj = 2eV/h being the Josephson angular frequency 
[16|), which allows us to use the quasi-static approach 
to treat the superconducting subsystem in contrast with 
the case analyzed in [lil [lj|. We demonstrate that a 
d.c. superconducting current may produce a strong ori- 
entation effect on the F layer magnetic moment. More 
interestingly, the a.c. Josephson effect, i.e. applying a 
d.c. voltage V to the iy9o-junction, would produce cur- 
rent oscillations and consequently magnetic precession. 
This precession may be monitored by the appearance of 
higher harmonics in CPR as well as a d.c. component of 
the current. In particular regimes, a total reversal of the 
magnetization could be observed. In the case of strong 
coupling between magnetic and superconducting subsys- 
tems, complicated non-linear dynamic regimes emerge. 

To demonstrate the unusual properties of the cpo- 
junction, we consider the case of an easy-axis magnetic 
anisotropy of the F material (see Fig[T]). Both the easy 
axis and gradient of the asymmetric spin-orbit potential 
n are along the z-axis. Note that suitable candidates for 
the F interlayer may be MnSi or FeGe. In these systems, 
the lack of inversion center comes from the crystalline 
structure, but the origin of broken-inversion symmetry 
may also be extrinsic, like in a situation near the surface 
of a thin F film. In the following, we completely disregard 
the magnetic induction. Indeed the magnetic induction 
in the (xy) plane is negligible for the thin F layer consid- 
ered in this paper, whereas the demagnetization factor 
cancels the internal induction along the z-axis (N = 1). 
The coupling between F and S subsystems due to the or- 
bital effect has been studied in [13] and it appears to be 
very weak, and quadratic over magnetic moment M for 
the case when the flux of M through the F layer is small 
in comparison with flux quantum $o = h/2e. 

The superconducting part of the energy of a ipo- 
j unction is 



Es {(f, (fio) ^ Ej[l- cos {ip ~ tpo)] , 



(1) 



where Ej = ^qIc/2tt is the Josephson energy, Ic is the 
critical current and ipo is proportional to the My com- 
ponent of the magnetic moment (see FiglT]). Therefore, 
when the magnetic moment is oriented along the z-axis, 
we have the usual Josephson junction with po = 0. As- 
suming the ballistic limit we may estimate the charac- 
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FIG. 1: Geometry of the considered (^o-junction. The ferro- 
magnetic easy-axis is directed along the z-axis, which is also 
the direction n of the gradient of the spin-orbit potential. The 
magnetization component My is coupled with Josephson cur- 
rent through the phase shift term ipo oc n. (M A V^*), where 
^ is the superconducting order parameter (V'l' is along the 
X-axis in the system considered here). 



teristic Josephson energy as [i] ^qIc/S Tckj^ sin £/£ 
with £ = AhL/hvp, where S, L and h are the section, the 
length and the exchange field of the F layer, respectively. 
The phase shift is 



nVsO My 

— TT 

vp Mo 



(2) 



where the parameter Vso/vp characterizes the relative 
strength of the spin-orbit interaction [4]. Further on we 
assume that Vgo/vp ~ 0.1. If the temperature is well be- 
low the Curie temperature, Mq = ||M|| can be considered 
as a constant equal to the saturation magnetization of the 
F layer. The magnetic energy contribution is reduced to 
the anisotropy energy 



E 



M 



KV f 

Mo 



(3) 



where K is an anisotropy constant and V is the volume 
of the F layer. 

Naturally, we may expec^ that the most interesting 
situation corresponds to the case when the magnetic 
anisotropy energy does not exceed too much the Joseph- 
son energy. From the measurements 18[ on permal- 
loy with very weak anisotropy, we may estimate K ~ 
4.10~^K.A . On the other hand, typical value of L 
in S/F/S junction is L 10 nm and sini/t ~ 1. Then, 
the ratio of the Josephson over magnetic energy would 
be Ej/Em ^ 100 for ~ 10 K. Naturally, in the more 
realistic case of stronger anisotropy this ratio would be 
smaller but it is plausible to expect a great variety of 
regimes from Ej/Em < 1 to Ej/Em > 1- 

Let us now consider the case when a constant current 
I < Ic is applied to the (po-junction. The total energy is 
(see, e.g. il6[): 



and both the superconducting phase shift difference ip 
and the rotation of the magnetic moment My = Mq sin 9 
(where 9 is the angle between the z-axis and the direction 
of M) are determined from the energy minimum condi- 



tions dipEtot 



dipfjEtot = 0. It results in 



-F with r 



Ej Vso 

KV vf' 



(5) 



which signifies that a superconducting current provokes 
the rotation of the magnetic moment My in the (jjz) 
plane. Therefore, for small values of the rotation, 9 (/) 
dependence is linear. In principle, the parameter F can 
be larger than one. In that case, when the condition 
I/Ic > l/F is fulfilled, the magnetic moment will be ori- 
ented along the y-axis. Therefore, applying a d.c. super- 
conducting current switches the direction of the magneti- 
zation, whereas applying an a.c. current on a (^o-junction 
could generate the precession of the magnetic moment. 

We briefiy comment on the situation when the direc- 
tion of the gradient of the spin-orbit potential is perpen- 
dicular (along y) to the easy axis z. To consider this 
case we simply need to take (po — £ (vso/vp) cos9. The 
total energy ([4]) has two minima 9 = (0,7r), while apply- 
ing the current removes the degeneracy between them. 
However, the energy barrier exists for the switch from 
one minimum into another. This barrier may disappear 
if F > 1 and the current is large enough / > /c/F. In 
this regime the superconducting current would provoke 
the switching of the magnetization between one stable 
configuration 9^0 and another 9 — n. This corresponds 
to the transitions of the junction between +ipo and ~ipo 
states. The read-out of the state of the </3o-junction may 
be easily performed if it is a part of some SQUID-like 
circuit (the (/jq-J unction induces a shift of the diffraction 
pattern by ipo). 

In fact, the voltage-biased Josephson junction, and 
thus the a.c. Josephson effect provides an ideal tool 
to study magnetic dynamics in a (^o-junction. In such 
a case, the superconducting phase varies with time like 
tp (t) — ojjt [l9|. If hujj ^ Tc, one can use the static value 
for the energy of the junction (jj]) considering ip (t) as an 
external potential. The magnetization dynamics are de- 
scribed by the Landau-Lifshitz-Gilbert equation (LLG) 

M 



— .^MxH.^-f-fMx- 



where Hcs — —SF/VSM. is the effective magnetic field 
applied to the compound, 7 the gyromagnetic ratio, and 
a a phenomenological damping constant. The corre- 
sponding free energy F = Es + Em yields 
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ipl + Es{ip,ipo) + Em {(po) , (4) where r = £vso/vF. Introducing rrn 
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(7) 
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FIG. 2: Results of numerical analysis of the magnetic moment dynamics of the (y9o-junction. a) Reversal of from analytical 
expression Eq. (|15|l (dashed curve) and numerical one (normal curve). The other curves are related to the M trajectory: b) 
in strong damping case c) and d) in the strong coupling regime revealing incomplete and complete magnetic moment reversal, 
respectively. 



(ujp = jK/Mq is the frequency of the ferromagnetic res- 
onance) in LLG equation ([S]) leads to 

rhx ~ rriz (t) uriy (t) — Trriz (t) sin (ujt — rmy) 
rhy = -m^ (t) (t) , (8) 

rhz — Trrix (r) sin {ujt — rmy) 

where uj — LUj/ujp. The generalization of Eq.(I5]) for 
a 7^ is straightforward. One considers first the "weak 
coupling" regime F ^ 1 when the Josephson energy Ej 
is small in comparison with the magnetic energy Em- In 
this case, the magnetic moment precess around the z- 
axis. If the other components verify {mx,my) <C 1, then 
the equations ([5]) may be linearized, and the correspond- 
ing solutions are 

rux (t) = '^"^ '^/^ and my (t) = - ^^^"^ '^''^ . (9) 

Near the resonance uj up, the conditions of lineariza- 
tion are violated and it is necessary to take the damping 
into account. The precessing magnetic moment influences 
the current through the (/jQ-junction like 

/ Tr 1 

— = sincjji + — ^5 sin2cj.7t + (10) 

ic 2 cj^ — 1 

i.e., in addition to the first harmonic oscillations, the 
current reveals higher harmonics contributions. The am- 
plitude of the harmonics increases near the resonance and 



changes its sign when ojj = ujp- Thus, monitoring the 
second harmonic oscillations of the current would reveal 
the dynamics of the magnetic system. 

The damping plays an important role in the dynamics 
of the considered system. It results in a d.c. contribu- 
tion to the Josephson current. Indeed, the corresponding 
expression for my it) in the presence of damping becomes 

my it) — smwjiH cosujjt, (11) 

r r 

where 

Fr w ± 1 Tr a 

with ri± = (w ± 1)^ + a^. It thus exhibits a damped 
resonance as the Josephson frequency is tuned to the fer- 
romagnetic one w — > 1. Moreover, the damping leads 
to the appearance of out of phase oscillations of my (t) 
(term proportional to coso^ji in Eq. pTjl ). In the result 
the current 

/ (t) « Ic sin u)jt -f Ic ^ ^ sin 2ujjt-{' 

+ /c "' ^ cos2cjj^ + /o(a) (13) 
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acquires a d.c. component 



(14) 



This d.c. current in the presence of a constant voltage V 
appHed to the junction means a dissipative regime which 
can be easily detected. In some aspect, the peak of d.c. 
current near the resonance is reminiscent of the Shapiro 
steps effect in Josephson junctions under external r.f. 
fields. Note that the presence of the second harmonic 
in / (t) Eq. (fT5|) should also lead to half- integer Shapiro 
steps in (/?o-junctions [21 1. 

The limit of the "strong coupling" T > 1 (but r < 1) 
can also be treated analytically. In this case, my « and 
solutions of Eq. ([8]) yields 



rrix (t) = sin 
niz (t) — cos 



r , 

— (1 — COSLU jt) 



(1 — cosojjt) 



(15) 



which are the equations of the magnetization reversal, a 
complete reversal being induced by T/lo > 7r/2. Strictly 
speaking, these solutions are not exact oscillatory func- 
tions in the sense that ttiz (t) turns around the sphere 
center counterclockwise before reversing its rotation, and 
returns to the position niz {t = 0) = 1 clockwise, like a 
pendulum in a spherical potential (see Fig[21c). 

Finally, we have performed numerical studies of the 
non-linear LLG Eq.® for some choices of the parame- 
ters when the analytical approaches fail. To check the 
consistency of our numerical and analytical approaches, 
we present in Figl^la) the corresponding ruz (t) depen- 
dences for low-damping regimes. They clearly demon- 
strate the possibility of the magnetization reversal. In 
Figsl2)3-d), some trajectories of the magnetization vec- 
tors are presented for general coupling regimes. These re- 
sults demonstrate that the magnetic dynamics of S/F/S 
t/jQ-junction may be pretty complicated and strongly non- 
harmonic. 

If the (/jQ-junction is exposed to a microwave radia- 
tion at angular frequency Wi, the physics that emerge 
are very rich. First, in addition to the Shapiro steps 
at ujj — ncji, half-integer-steps will appear. Secondly, 
the microwave magnetic field may also generate an addi- 
tional magnetic precession with ui frequency. Depend- 
ing on the parameters of (/?o-junction and the amplitude 
of the microwave radiation the main precession mecha- 
nism may be related either to the Josephson current or 
the microwave radiation. In the last case the magnetic 
spin-orbit coupling may substantially contribute to the 
amplitude of the Shapiro steps. Therefore, we could ex- 
pect a dramatic increase of this amplitude at frequencies 
near the ferromagnetic resonance. When the influence 
of the microwave radiation and Josephson current on the 



precession is comparable, a very complicated regime may 
be observed. 

In the present work we considered the case of the easy- 
axis magnetic anisotropy. If the ferromagnet presents an 
easy-plane anisotropy than qualitatively the main conclu- 
sions of this article remain the same because the coupling 
between magnetism and superconductivity depends only 
on the My component. However, the detailed dynamics 
would be strongly affected by a weak in-plane anisotropy. 

To summarize, we have demonstrated that S/F/S ipo- 
junctions provide the possibility to generate magnetic 
moment precession via Josephson current. In the regime 
of strong coupling between magnetization and current, 
magnetic reversal may also occur. These effects have 
been studied analytically and numerically. We believe 
that the discussed properties of the iy9o-junctions could 
open interesting perspectives for its applications in spin- 
tronics. 

The authors are grateful to Z. Nussinov, J. Cayssol, 
M. Houzet, D. Gusakova, M. Roche and D. Braithwaite 
for useful discussions and comments. This work was sup- 
ported by the French ANR Grant N° ANR-07-NANO- 
011: ELEC-EPR. 



[4] 
[5] 
[6] 
[7] 
[8] 

[9] 
[10] 

[11] 
[12] 
[13] 
[14] 

[15] 

[16] 

[17] 
[18] 
[19] 



Electronic address: |f.konschelle@cpmoh.u-bordeauxl.fT| 

Also at Institut Universitaire de France. 

I. Zutic, J. Fabian, and S. Das Sarma. Rev. Mod. Phys. 

76, 323 (2004). 

J. Hauptmann, J. Paaske, and P. Lindelof. Nature 
Physics 4, 373 (2008). 

R. Winkler. Spin-orbit coupling effects in two- 
dimensional electron and hole systems (Springer, New 
York, 2003). 

A. Buzdin. Phys. Rev. Lett. 101, 107005 (2008). 

A. A. Reynoso et al. Phys. Rev. Lett. 101, 107001 (2008). 

A. L Buzdin. Rev. Mod. Phys. 77, 935 (2005). 

C. Bell et al. Phys. Rev. Lett. 100, 047002 (2008). 

J.-X. Zhu, and A. V. Balatsky Phys. Rev. B 67, 174505 

(2003). 

L. Bulaevskii et al. Phys. Rev. Lett. 92, 177001 (2004). 

J.-X. Zhu et al. Phys. Rev. Lett. 92, 107001 (2004). 

Z. Nussinov et al. Phys. Rev. B 71, 214520 (2005). 

S. Takahashi et al. Phys. Rev. Lett. 99, 057003 (2007). 

M. Houzet. Phys. Rev. Lett. 101, 057009 (2008). 

X. Waintal, and P. W. Brouwer. Phys. Rev. B 65, 054407 

(2002). 

V. Braude, and Y. M. Blanter. Phys. Rev. Lett. 100, 
207001 (2008). 

K. K. Likharev. Dynamics of Josephson junctions and 
ctrcmts (Gordon and Beach Science Publishers, 1986). 
S. Hikino et al. J. Phys. Soc. Jap. 77, 053707 (2008). 

A. Y. Rusanov et al. Phys. Rev. Lett. 93, 057002 (2004). 

B. D. Josephson. Superconductivity (in two volumes), 
vol. 1, chapter 9, Weakly coupled superconductors (R.D 
Parks, Marcel Dekker, Inc., 1968). 



E. M. Lifshitz, and L. P. Pitacvskii. Course of theoret- [21] H. Sellier et al. Phys. Rev. Lett. 92, 257005 (2004). 
teal physics, Volume 9 : Theory of the condensed state 
(Butterworth Heinemann, 1991). 



